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We argue that the strong jet quenching of heavy flavors observed in heavy-ion coUisions is to a 
large extent due to binary scatterings in the quark-gluon plasma. It can be understood from first 
principles: the charm collision probability beyond logarithmic accuracy and Markov evolution. 

PACS numbers: 12.38Mh 
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In particle physics, heavy flavors are suitable to inves- 
tigate properties of the strong interaction. Their large 
mass also proves useful when studying aspects of many- 
body QCD for heavy-ion phenomenology, since quarks 
with TO 3> T can be considered as test particles in a 
heat bath. Given this simplification, heavy quarks (say 
charm with TOc « 1.2 GeV, which might be sufficiently 
large for temperatures reached at RHIC) can help to clar- 
ify the ongoing debate on the parton energy loss in the 
quark-gluon plasma. For light partons, it is common to 
attribute jet quenching entirely to radiative energy loss, 
motivated by its parametric energy dependence [1, 2]. 
One may note that plasma parameters inferred in such 
approaches can be hard to reconcile with general expec- 
tations, cf. [3]. What is more, the pureJy radiative picture 
of jet suppression has recently been challenged by the ob- 
servation that heavy quarks (which radiate less [4]) are 
quenched almost as much as light quarks, as concluded by 
analyzing electron yields from heavy flavor decays [5, 6]. 
This has revived an interest [7, 8, 9, 10] for the collisional 
energy loss as an additional suppression mechanism in 
the range of moderately large momenta. 

With regard to the importance of the understanding 
of jet quenching, it seems necessary to point out that ex- 
isting approaches suffer from several fundamental short- 
comings, which makes their conclusions evasive. A ma- 
jority of approaches is based on Fokker-Planck equations, 
despite the fact that the assumed dominance of soft scat- 
tering [11] is justified only at leading logarithmic accu- 
racy [12]. As one consequence, detailed balance must 
be imposed in some way, whereas equilibration should 
be a prediction of the formalism. A second severe is- 
sue of existing approaches is the disregard of the mo- 
mentum dependence of the strong interaction. Strictly 
speaking, QCD calculations are not predictive without 
discussing the running coupling (opposed to QED with 
Q!qjj,j^ ~ 1/137 in the widely applicable Thomson limit). 
For quantitative estimates, a value of the coupling has 
then to be assumed; common is afix ~ a(— (27rT)^). 
However, the energy loss probes also parametrically soft 
scales, ~ -v/a T, which is not only decisive for its high- 
energy limit [22]; it will also turn out to be crucial for 
heavy ion phenomenology. 

Let us start by considering the average energy loss of a 
charm quark subject to binary collisions with gluons and 



light quarks (i = g, q) in the thermalized plasma, 
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Here d = 6 for a charm quark; v = p/ E is its velocity, 
ni{k) = Xjie^l^ =F 1) and = 1 ± Jii are the distribution 
functions of the collision partners, and lo = E — E' \& the 
energy transfer in the scattering. \M.i\ is summed over 
color and spin states of all particles. 

Evaluating Eq. (1) with the Born cross sections [13] 
would yield divergent results because of long-range gauge 
interactions. Therefore, screening effects (which formally 
arise from thermal loop corrections) have to be taken into 
account already in a leading order calculation. These 
thermal corrections come along with the vacuum fluc- 
tuations, which diverge and need to be renormalized. 
Renormalizability, as such, is obviously not affected by 
the (UV-finite) thermal contributions - which might be a 
reason why renormalization is often utterly disregarded 
in finite-temperature field theory. However, only by this 
fundamental concept QCD calculations can be predictive. 

Rigorous renormalization can be tedious, especially in 
thermal field theory. Fortunately, for the observables of 
interest here, the leading-order results can be inferred 
by elementary reasoning. For gluon exchange processes 
(t-channel scattering) it has been argued previously [14] 
that resumming and renormalizing the loop corrections 
amounts to replacing, in the Born amplitudes, the bare 
coupling by the running coupling a{t)^ schematically 
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where q = (w^ — <)^/^. We will parameterize the thermal 
self-energy Ht by an effective cut-off, of the order of the 
Debye mass and evaluated with running coupling [15], 



p'^{t) = K • 47r (l 



a{t) 
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The customary ad hoc choice for this screening mass is 
/^ix = ^TT (l -I- i n/) Q!fix T^, which we adopt only for the 
sake of comparison with existing estimates. In fact, the 
coefficient k in Eq. (3) can be calculated by compar- 
ing c?i5*~'^'^™"°Yc?a;, evaluated with effective cut-off, to 
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the strict result beyond logarithmic accuracy [16], which 
yields 

K = (2e)"i « 0.2. (4) 

The resulting cut-ofF is small, which will be a main source 
of differences between our and prevalent estimates. 

Complementing the arguments put forward in [14] , we 
now turn to the s and u-channel contributions to charm 
scattering. In order to quantify 'the' coupling, consider 
the vacuum corrections to the Born amplitudes, which 
have logarithms of characteristic momenta. By specify- 
ing a renormalization scale /Uij, bare quantities are ex- 
pressed in terms of physical ones, such as the coupling 
a(/i|j.). By renormalization flow equations, all choices 
of fifj are physically equivalent. This allows, if there is 
only one characteristic momentum P, to choose = P. 
Then the logarithmic term vanishes, and the renormal- 
ized result looks like the Born approximation - albeit 
with a a{P^). Now, the coUisional energy loss (1), as 
a thermal average of |A^i|^, is dominated by interactions 
with 

this particular kinematics, there is indeed only one rele- 
vant momentum scale determining a - the virtuality 
of the intermediate state, cf. [17]. We will fix the couphng 
also for the remaining sub- dominant scattering contribu- 
tions at the respective virtuality, not without estimating 
the arising uncertainty. For time- like contributions, we 
continue the 1-loop coupling according to Ref. [18], 
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FIG. 1: Probability density E,=g,9 -P»('^.P = 8 GeV). The 
band gives the sensitivity under variation u G [^,2] of the 
scale in the running coupling a{vV'^). Shown for comparison 
is the result obtained with the coupling qax and cut-off /iflx. 

for a test particle with momentum p to change, per unit 
of time, its energy by uj. The condition H > 0, where 
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for g2 ^ , (5) 



where /3o = H^f with nj = 3, and L± = ln(±Q^/A^). 
With regard to quantitative estimates it is worthwhile 
recalling that perturbative approaches can be of use at 
surprisingly soft momentum scales [19]. May details of 
the behavior of a{Q^) in the deep infrared be uncertain, 
there exists a robust constraint (universality hypothesis) , 



Aq^\<qi 
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where Qu = 2 GeV [19]. Accordingly, we impose an up- 
per bound on the running coupling (5), a{Q^) < 1.1 for 
our preferred QCD parameter A = 0.2 GeV (adjusted to 
lattice results for the heavy quark potential, cf. [15]). We 
have verified that the deep infrared region as well as the 
precise value of A are not very important for our concerns 
~ for s and u-channel terms due to phase space suppres- 
sion, and for ^-channel contributions by screening. 

Returning to Eq. (1), the thermal 2-body phase space 
can be expressed as an integral over lo and the invariant 
momentum transfer t [12]. This allows us to write 

^ dui Pi{uj,p)uj , 



determines not only the bounds t± in Eq. (7), it also 
constrains the integral over k, in particular such that 
fc -|- cj > is always fulfilled. Thus energy gain is expo- 
nentially suppressed for uj ^ —T. For lo > P ^ lo^'^ 
up to near the kinematic threshold, see Fig. I , which also 
shows the uncertainties arising from the scale setting. 

In the following we consider an ensemble of test par- 
ticles. In order to describe its momentum distribution 
fij)) it is convenient to introduce the rate ^^{p'jp) for 
transitions p ^ p' . The spectrum changes by scattering 
out and into a given state: for a sufficiently small time 
interval f{q,St) = {1-Str{q)) f{q) + 6tJdpV{q,p) f{p), 
where r{q) ~ J dp'P{p, q) denotes the interaction rate of 
particles with momentum q. Introducing 

T{q,p) = (1 - StT{q)) 5{q-p) + StV{q,p) (9) 

makes explicit that the time evolution of the spectrum. 



f{q,t + St)^ dpTiq,p)f{p,t), 



(10) 



is a first order Markov process (the transition depends 
only on the state rather than on preceding history). We 
discuss the resulting Markov chain on a discrete momen- 
tum space, which arises naturally from binning [24] and 
reduces the convolution ( Id) to a matrix multiplication. 
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dx 

with the probability density [23] 
ni{k) 



P.ito,p) 



2k 



■ni{k + Lu) 
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In other words, the evolution in (discrete) time is de- 
termined simply by powers of the transition matrix Tqp, 
which is not only conducive to numerical studies, it also 
brings instructive insight readily [25]. Indeed, essential 
properties of the evolution follow alone from 
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which holds by definition of r{q). It imphcs first that Tqp 
maps the hyperplane H characterized by ||/|| = fp = 
constant onto itself. Put differently, the 'norm' ||/|| is in- 
variant under Tqp, which in the present context describes 
particle number conservation. From this we infer |Ai| < 1 
for the eigenvalues of Tqp-, otherwise the norm would not 
be conserved in repeated mappings. Among {A^} is an 
eigenvalue Ai = 1 since the rank of the matrix Tqp — 1 
is dim(7^j,) — 1. These properties of the eigenvalue spec- 
trum arc proven more rigorously in the Perron- Frobenius 
theorem, cf. [20], which also shows that Ai is a simple 
eigenvalue. Consequently, the stationary state 

rq^^ hm r;;/p(o), (13) 

is unique and approached for each initial distribution. 

The corresponding eigenvectors of Tqp form a basis 
with the peculiarity that only ei has a component per- 
pendicular to the hyperplane, i.e., e J? for i ^ I, 
see Fig. 2. Hence, we can separate in the evolution 
fin St) = E.=iA?</>.e, = + E,=2 A? 0»e, a part 
which proceeds entirely in H. Since |Ai^i| < 1, the ap- 




FIG. 2: Equilibration in 2 dimensions (where H is given by 
x + y — 1); the initial states are (1,0) and (0,1), respectively. 

proach to equilibrium for large t is exponentially fast, 

\\fit)-n\\ (X exp(-t/r--), (14) 

where the so-called mixing time 

r'"''^ = (5i/ln(A^i) (15) 

is determined by the second largest eigenvalue [26]. 

It is instructive to illustrate these essential features of 
thermalization by means of a toy model, namely 

Tij ^ {1- I - g) Stj + I (5ij+i - g (5ij_i , (16) 

with constants I and g parameterizing the relevant loss 
and gain rates in Tqp. First, it is easy to see that the 
stationary distribution is exponential, f^^ oa exp(— 1/0), 
with the 'temperature' Q = l/hi(l/g) being determined 
by the logarithmic ratio of loss and gain [27] . Further- 
more, for large dimension of Tj, one can readily derive 

\2 = l-iVl-^)\ (17) 

which shows a non-analytic behavior of the mixing time. 

In QCD, for sufficiently fast charm quarks, the 
(binned) loss rate dominates over gain (cf. Fig. »), thus 

r""'' ~ T^^l^ - F^;;,^^! . (18) 



This relation between long and short-time aspects of 
equilibration is worth emphasizing. More important. 
Fig. 3 reveals that conventional calculations would over- 
estimate this crucial time scale by a factor /-C « 5 ~ 
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FIG. 3: Average time t^^ for a charm quark to change (gain 
or lose) momentum by Ap > 0.4 GeV, compared to the fixed- 
coupling estimate (uncertainty band as in Fig. 1). 

Realizing that binary collisions are far more effective 
than previously estimated, we consider in an exploratory 
study jet quenching in the mid-rapidity region of central 
collisions within the Bjorken model. To compare with ex- 
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TABLE I: Representative parameter sets for the Bjorken 
model (units: GeV or fm; *adjusted to total entropy S ~ 10"*). 



isting results [8, 9], we adopt parameterizations as sum- 
marized in Tab. I. For transparency, we will not discuss 
effects of hadronization (for heavy flavors the partonic 
suppression can be indicative for the observed Raa)- 

Consider first the evolution of fo{Pt) oc ptdNini/dpf 
at constant temperature. Fig. 4 shows the hard part of 
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FIG. 4: Evolution of the initial charm spectrum / at fixed T 
(for uncertainty band cf. Fig. J). 

the spectrum being quenched markedly already after a 
few fm/c, as could be anticipated from r^p ~ Ifm/c. 
Taking then into account the path length distribution 
dN/dl cx (1 - (Z/2i?)2)i/2 in Bjorken's heat bath, T = 
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Tq{t /tq)~^/^ , is straightforward in the Markov formal- 
ism. Fig. 5 shows the eharm quenching ratio taa ~ 
{dN{t\{ie)/dp1)/{dN\ni/dpl)\ to compare to taa ^ 0.85 
with fixed couphng, and to Refs. [7, 8, 9, 10]. Given that 
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FIG. 5: Comparison of partonic charm quenching in Bjorken 
model to the observed Raa- For better visibility, the inter- 
section of uncertainty bands (cf. Fig. I ) was not hatched. 

we have calculated the underlying mechanism from first 
principles, the comparison with the RHIC data [5, 6] is 
striking. Depending on details of the initial spectrum 
and hadronization, there is room for radiative quench- 
ing (with reasonable parameters), which becomes more 
important at larger pf With regard to the increasing 
Taa at smaller pt we underline that the Cronin effect 
has not been taken into account in the calculation. The 
small value of rAA at ^ 4 GeV is quite robust, as veri- 
fied here with two rather different parameter sets within 



the Bjorken model. Radial expansion will reduce the life 
time of the plasma phase, but also modify the path length 
distribution dN/dl. The resulting compensation will be 
quantified in a forthcoming study with more realistic col- 
lision dynamics and including hadronization [21]. 

To summarize, we have shown in a Markov formalism 
that, for moderately large momenta, binary collisions are 
a key mechanism for jet quenching of heavy flavors. This 
conclusion requires taking into account - more carefully 
than in existing approaches - essential features of QCD, 
namely the momentum dependence of the strong cou- 
pling and relevant screening effects. Both aspects are 
interconnected and, as a matter of fact, mandatory for 
quantitative estimates from thermal field theory, unless 
temperatures are asymptotically large. 

In closing, since jet quenching has significantly influ- 
enced the current picture of the strongly coupled quark- 
gluon plasma (sQGP), it seems worthwhile taking a more 
general point of view. Commonly used but rather crude 
estimates (as for the relevant screening range) can lead to 
interpretations of a 'too' strongly coupled plasma. While 
for heavy ion phenomenology the relevant coupling is cer- 
tainly not small, we do see the possibility to understand 
essential features of the quark-gluon plasma from flrst 
principles - perturbative QCD can work like a charm. 
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For i? — > oo, the energy loss from Eq. (1) is of first order 
in the coupling [14], notwithstanding \Mi\^ ^ o? ■ 
Note that only the first moment of PiuS) is logarithmi- 
cally enhanced, which justifies the parameterization (.i). 
Binning regulates the divergence of the interaction rate 
r(g) = J dpVip, q) arising from p ~ g. 
It will give a supplementary perspective on conclusions 
derivable from Boltzmann's equation. 
This consideration gives a transparent interpretation of 
the relaxation time in the homonymous approximation. 
The charm transfer matrix Tqp will of course lead to equi- 
libration at the temperature of the heat bath, see Fig. 4. 



